<Z3 



Practical Methods for Wireless Network Coding 
with Multiple Unicast Transmissions 

Tugcan Akta§, A. Ozgiir Yilmaz, Emre Akta§ 



(N 

q . Abstract 

We propose a simple yet effective wireless network coding and decoding technique for a multiple unicast 
network. It utilizes spatial diversity through cooperation between nodes which carry out distributed encoding 
operations dictated by generator matrices of linear block codes. In order to exemplify the technique, we make 
use of greedy codes over the binary field and show that the arbitrary diversity orders can be flexibly assigned to 
nodes. Furthermore, we present the optimal detection rule for the given model that accounts for intermediate node 
errors and suggest a low-complexity network decoder using the sum-product (SP) algorithm. The proposed SP 



. detector exhibits near optimal performance. We also show asymptotic superiority of network coding over a method 

> 

that utilizes the wireless channel in a repetitive manner without network coding (NC) and give related rate-diversity 

o" 

. trade-off curves. Finally, we extend the given encoding method through selective encoding in order to obtain extra 



coding gains. 
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I. Introduction 

In order to counteract the effects of fading in wireless communication networks, many ways of creating 
diversity for transmitted data have been proposed. Utilizing the spatial diversity inherent in wireless 
channels, cooperative communication [Q]| has been of great interest in recent years. In [0, [0 three 
methods to be used by relay nodes are described: amplify-forward (AF), decode-forward (DF) and detect- 
forward (DetF). The AF method attains full diversity, whereas other two cannot, unless the propagation 
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of errors resulting from the decoding/detection operation is avoided, for example by using a selective 
transmission strategy that advocates the forwarding of only reliable enough packets. One of the various 
ways to handle this problem is using CRC -based methods, which results in loss of spectral efficiency due 
to drop of a packet with only a few bit errors. An on/off weighting based on relay signal-to-noise power 
ratio (SNR) is given in 01. Weighting of the signals either at the relay or at the receiver using the relay 
error probability is proposed in Q, 10. Yet another idea is transmitting the log-likelihood ratios (LLR) 
of bits 0. However, the soft information relaying methods in 0-0 suffer from quantization errors 
and high peak to average ratio problems. In addition, the AF method requires hardware modifications on 
modern-day communications systems and the DF method leads to high complexity decoding operations 
especially for the relays. As an alternative, relays may use the simple DetF method, which is shown to 
avoid error propagation in [3 J, if the error probabilities at relays are known and the maximum a posteriori 
probability (MAP) detection is employed at the receiver. 

NC was initially proposed to enhance network throughput in wired systems with error-free links of unit 
capacity [8]. Later studies exhibited the good performance of random linear NC [9]. In wireless networks 
with nodes naturally overhearing transmissions, NC can be utilized to create diversity, reduce routing 
overhead, and introduce MAC layer gains as discussed for practical systems in IfTOll . Although most of 
the work in the literature concentrate on multicast transmission ifTTI : we deal with a network involving 
multiple unicast transmissions, which is inherent in real-life scenarios. Hence we formulate a multiple 
unicast transmission problem such that for each unicast transmission, there is a distinct diversity that is 
to be improved via spatial opportunities. 

We consider a simple NC scheme based on DetF. Given a relay combining strategy, which we represent 
by a generator matrix and a vector of transmit schedule, we investigate the diversity order of each source, 
which can be unequal. We propose a novel method for designing the generator matrix based on linear 
block codes over the binary field. The proposed method is very flexible in that any set of desired diversity 
levels for the sources can be achieved with the highest NC rate possible. The diversity analysis relies on 
an optimal MAP decoder at the destination which employs the reliability information of the intermediate 
nodes and avoids loss of diversity due to error propagation J3]|, |fl2ll . The numerical complexity of the given 
optimal decoder can be impractical. Thus we propose a practical approximation of the MAP detector: the 
SP network decoder. 

A study based on flexible network codes in a two-source two-relay system with emphasis on unequal 
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error protection is |fT3lh where authors propose a suboptimal detection rule (distributed minimum distance 
detector) that is known to result in diversity order loss. Our scheme captures full diversity due to the 
use of the SP detector with intermediate node reliability information. In [14], a multicast scenario is 
investigated (again under additive white Gaussian noise (AWGN) assumption) for obtaining an optimal 
energy allocation scheme in order to minimize bit error rate at the sink nodes. In [fT5l . performance of a 
multiple hop network with no fading assumption is analyzed in terms of effects of the bit errors at the 
relays utilizing a technique known as error event enumeration. Similarly in [fT6ll , optimal detection rule 
to be used at the destination node is presented for AWGN channels in addition to the description of a 
genie-aided decoder which yields a lower bound on the performance of the optimal detector. Different 
than |fT5l and ffTBTh we consider faded wireless links and give optimal detection rule corresponding to 
this realistic scenario. Futhermore, the operation at the intermediate nodes in our scheme is DetF as 
opposed to the more complex DF in ffT5lh [fT6l . One of the studies closest to ours is IfPTIh where the 
NC operation is fixed in construction yielding very large Galois field (GF) sizes for increasing network 
size and relay nodes carry out complicated DF operation for each transmission they overhear. Similarly 
in |[T8lh DF is used in a fixed single-relay two-user scenario in order to provide diversity-multiplexing 
trade-off for NC. However, our results indicate that any diversity order can be achieved for any unicast 
transmission even with the GF of size 2 by using linear block codes as the network codes and simple DetF. 
Also independent from our work, in lfT2l . [fP9l , results concerning diversity analysis for a system model 
resembling ours have been obtained. Similarly in a recent work ll20l . additional coding gain analysis is 
given for a multi-source multi-relay network with relays having no data to be transmitted. Our model is 
generalized in the sense that each node both acts as a source node with its own data to be conveyed over 
the network and as an intermediate node serving as a means for combining and relaying others' data. 
Moreover, the proposed detection rules of lfT2ll . [fP9ll , [|20l result in exponential decoding complexity in 
the number of transmissions, since they are based on maximum likelihood sequence estimation. Recently 
in ETTl . a wireless broadcast network with block erasures is considered and a network coding scheme 
is proposed for retransmissions. The improvement in the number of retransmissions for the downlink 
channel with respect to the conventional automatic-repeat-request mechanism is clearly shown. Another 
recent work [22] identifies the diversity-multiplexing trade-off for a NC system, in which multiple access 
to the channel is allowed to be non-orthogonal. On the other hand, our model with orthogonal access of 
the nodes to the channel does not require a complex successive interference cancellation technique to be 
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implemented at the destination node. 

The major goal in this paper is to introduce practical NC/decoding methods for improving the diversity 
order of a network through cooperation with the overall rate of transmission in mind. The contributions of 
the current paper can be listed as follows. (1) A generalized wireless NC scenario with nodes possessing 
both relay and source roles and its diversity analysis. (2) Design of novel network codes based on close-to- 
optimal linear block codes. (3) Investigation of greedy codes and maximum code rates for desired diversity 
levels. (4) Application of the SP algorithm for decoding network codes with relay reliability information. 
The rest of the paper is organized as follows. In Section [III we present the wireless network model, the 
corresponding detection rule that is individually optimal for each user, and a practical enhancement on 
the proposed network encoding method. We investigate the details of network codes based on linear block 
codes with emphasis on greedy NC through some sample networks in Section [Till. We also give some 
asymptotic results based on the rate-diversity order trade-off curves for the proposed NC method and a 
repetitive method that represents the no NC case. In Section [TV] we introduce the SP network decoder 
that has linear complexity order and yields performance figures very close to that of the optimal decoder. 
Section [V] includes the numerical results for the mentioned network encoding and decoding methods and 
Section [VTj concludes the paper. 

II. Wireless Network Model 
A. The Network-Coded System and The Corresponding Separation Vector 

In this work, we analyze a wireless network in which unicast transmission of data symbols, each 
belonging to a different source, is to be carried out utilizing NC at the intermediate nodes. Under the 
general operation scheme, every node may act both as a member (source or destination) of a unicast 
communication pair and as an intermediator (relay) node for other unicast pairs. 

In order to ease the explanation of system model and the roles of nodes in the network, here we 
start a with a simple network-coded operation depicted in Fig. [TJ which makes use of binary NC 
through usual binary addition operator ©. The network of interest consists of k = 3 source nodes 
and a dummy node that represents a hypothetical detector of source packets u 1; u 2 , and u 3 at the 
corresponding destination nodes. The transmission of the these 3 data packets is allowed to be completed 
within n = 4 orthogonal time slots, which form a round of NC communication with a data rate of 
r = - = | packets /transmission slot. The channel is assumed to be shared by a time division multiple 
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access technique for the sake of simplicity in model description and due to causality requirements forcing 
the intermediate nodes to listen to a symbol before combining it through NC. 

As seen in Fig. [Q the first time slot is reserved for node 1 to transmit its own data packet ui and this 
transmission is overheard by source nodes 2 and 3 in addition to the destination node 0. We assume that 
the links between different pairs of nodes are independently Rayleigh faded. The channel corresponding to 
a link is quasi-static, i.e., constant over a packet and independently faded for different packets. We further 
assume that there is no feedback of channel state information (CSI) within the system in order to simplify 
the implementation and that each receiving node, including node 0, has the perfect knowledge of only 
the incoming links through measurements of the respective channels. Hence following the transmission 
of the packet ui, both node 2 and node 3 use respective CSI to obtain their own detection results on the 
symbols transmitted. Due to the block fading assumption, we consider a single data symbol u% and its 
detection/decoding event representing all symbols in the packet. Therefore, corresponding to the detections 
of Mi, each node has also a reliability information based on the probability of error in the detection, which 
is in fact only a function of its own channel measurement result. In order to counteract the effects of error 
propagation, this reliability information is passed to the destination node by an intermediate node (node 
2 or 3), whenever it combines the detected symbol shown by u\ with its own and other nodes' symbols. 

In the second time slot, node 2 transmits its own symbol and this transmission is observed by all other 
nodes as well. In the following slot, a NC operation is carried out by node 3, which simply combines its 
own symbol and its detection result for the first time slot U\. In order to inform the destination node 0, 
node 3 has to append the error probability for the network coded symbol u\ © to the packet it formed. 
In the last slot, once again node 2 uses the channel to transmit the network encoded data «i © u 2 with 
its own estimate of U\ and appends the corresponding reliability information to the transmitted packet. 
Hence the destination node knows only the reliability information for the last two transmissions which 
incorporate NC, but not the CSI between all intermediate nodes within the system. For the proposed 
scenario, the overhead of appending reliability information to the network-coded packets on the spectral 
efficiency is small for large packet lengths. Therefore, the effect of sharing the the reliability information 
on spectral efficiency is expected to be small. 

Up to this point, the sample NC system is detailed in words. From the perspective of destination node 
0, the same system can be described formally using a generator matrix G (called the transfer matrix 
in ifTTIO and a scheduling vector v. The columns of G, g'-s, represent the combining operations at the 
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intermediate nodes and the entries of v give the scheduling of the nodes accessing the channel: 



G 



10 11 
10 1 
10 



v=[1 2 3 2]. (1) 



The combined data vector u = [u\ w 2 u 3 ] can then be used to form the distributed codeword: c = 
[ci c 2 c 3 c 4 ] = uG. The choices u, G, k, and n for the parameters defining the operation of network are 
not arbitrary. They are used intentionally to point out the analogy to regular linear block codes. However, 
reliable detection of all data symbols, i.e., whole block u, originating from a single error-free source is 
of interest for a regular decoder; whereas node may desire to reliably detect, as an example, only ui 
under cooperative encoding. Hence we need to identify a parameter that describes the performance for 
detection of a single symbol u\ as opposed to the codeword u for our model. 

One can show that the minimum distance for G is 1 . However, we will see that an error event requires at 
least 2 bit errors for decoding of U\ at node 0. Let all the data bits be equal to without loss of generality, 
i.e., u = [0 0]. Hence the transmitted codeword is expected to be c = [0 0] for the case of no 
intermediate node errors. The error event for U\ corresponds to its decoding as 1. This erroneous decoding 
can occur for sequence decoding u G {[100], [101], [110], [111]}, where u denotes the decoding result at 
node 0. The incorrect codewords c corresponding to these decoded vectors are [1011], [1001], [1110], [1100], 
respectively. When these codewords are compared to the codeword [0000], it is clear that at least 2 bit 
errors are needed to cause an error event. Hence the minimum distance for u\ in this setting is said 
to be 2. The erroneous decoding for other bits can be investigated in a similar fashion. Focusing on 
ii 3 and hypothesizing u = [0 0], w 3 is incorrectly decoded when u G {[001], [011], [101], [111]}- The 
corresponding codewords are [0010], [0111], [1001], [1100]. Therefore, a single bit error can cause erroneous 
decoding of u 3 yielding a minimum distance of 1. As seen in the example, the error performance varies 
from symbol to symbol. Next we generalize this claim to cover arbitrary generator matrices and verify it 
through simulations in Section IVl 

Now we consider a subset of nodes in which there are k nodes transmitting data to a single destination 
node 0. Let the symbol transmitted by node i be denoted by for i G {1, ...,&}, and m be an 
element from the Galois field of size q, GF(g). We assume m to be statistically independent and define 
u = [ui u 2 ... ttfc] as the combined data vector. In time slot j G {l,...,n}, a transmitting node 
Vj G {1, . . . , k} forms a linear combination of its own and other nodes' data. If Vj detects all data to 
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be encoded correctly, it simply forms Cj = ugj, where gj is a k x 1 network encoding vector whose 
entries are elements of GF(g). Let u,i denote the estimate of the symbol of node i at node Vj. Using these 
estimates, node Vj forms the noisy network coded symbol dj = ugj that is also an element of GF(g). 
Then Vj modulates and transmits this symbol to receiver node as: 

Sj = fi(cj), (2) 

where shows the mapping of a symbol to a constellation point. Although symbols may come from 
any alphabet and non-binary constellations may be used, we will focus hereafter on GF(2) and binary 
phase-shift keying (BPSK) with Sj = 1 — 2c j. Our assumption is that vector gj, source address Vj and 
probability of error p e . for the transmitted symbol are appended to the corresponding packet. We consider 
transmissions with no channel coding and deal with single network coded data symbol Cj which represents 
all symbols within a packet transmitted by vj. At the end of a round of transmissions, if no errors occur 
at the intermediate nodes, the overall vector of n symbols coded cooperatively in the network is 

C = [d C 2 ... C n ] = U [gi g 2 ... g n ] = UG. (3) 

The generator matrix characterizes the network code together with the vector of transmitting nodes 

v = [vi v 2 ... v n ) . (4) 

Equations © and © generalize the definition of the example NC in (OQ). Next, we present Algo- 
rithm [Q which generalizes the method for finding the minimum distance for U{. In Algorithm \T\ the 



Algorithm 1 Algorithm for finding the minimum distance corresponding to symbol Ui for a (n, k, d) code 
with given generator matrix G 

minimumdistance n 
indexvector <— [1 2 ■ ■ ■ k]\i 
for j = 1 to 2 k - 1 do 

errorpattern <— dec2GF2(j — 1) 

errdatavector[i] ^— 1 

errdatavector [indexvector] errorpattern 
errcodevector <— errdatavector * G 
err codedi stance number ofnonzero(errcodevector) 
minimumdistance <— min(minimumdistance, errcodedistance) 
end for 



function dec2GF2(.) returns a binary pattern corresponding to the input decimal number and the function 
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number ofnonzero(.) returns the number of non-zero entries in the input vector. It is assumed that the 
data vector u consists of all Os, relying on the linearity of the network code. The algorithm first creates 
all possible erroneous data vectors u that have 1 in the ith position so that all possible codewords leading 
to erroneous decoding of are generated by c = uG. Afterwards, we search within these codewords 
to find the one with the minimum distance to the transmitted codeword of all Os. This minimum value 
gives us the minimum distance for U{. The set of minimum distances corresponding to all u^s (named 
the separation vector 11231 ) is utilized in identifying the performance metrics for NC in Section IIII-AI 

B. Optimal Network Decoding Using Intermediate Node Reliability Information 

The intermediate nodes are assumed to use the DetF technique (hard decision with no decoding 
operation) due to its simplicity. In a wireless network, an intermediate node Vj has a noisy detection 
result u of u. Let us express the resulting noisy network coded symbol as 

Cj — Cj © 6j, (5) 

where e, denotes this propagated error. We observe that a possible error in u propagates to dj after the 
NC operation dictated by gj is realized. We assume that node Vj knows the probability mass function 
(PMF) of 6j, p(ej), which we name as the intermediate node reliability information. This assumption is 
not unrealistic as it can be determined by the estimation of the channel gains of the links connected to 
Vj, along with the reliability information forwarded to Vj. The received signal by node at time slot j is 
then yj = hjSj + Wj, where hj is the channel gain coefficient resulting from fading during the jth slot 
and Wj is the noise term for the link between Vj and node 0. The gain coefficient is circularly symmetric 
complex Gaussian (CSCG), zero-mean with variance E s , i.e., it has distribution CN(0,E s ). The noise 
term is CSCG with CN(0, N ). The usual independence relations between related variables representing 
fading and noise terms exist. The overall observation vector of length n at node is 

y = Hs + w, (6) 

where y = [y 1 ... y n ] T , s=[s\ ... s n ] T = /i(c T ), w = [w% ... w n ] T and H is a diagonal matrix whose 
elements are independent channel gains hi, h 2 , . . . , h n for the links connected to node 0. It is assumed 
that H is perfectly known at node 0. Combining the coded symbols in a network code vector, we obtain 



c = c©e = uG©e, 



(7) 
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where e = [ei ... e n ] is the error vector denoting the first hop errors. We assume that e is independent of 
c although dependence can be incorporated in the SP decoder developed in Section [IV] This independence 
assumption is valid directly for BPSK modulation, whereas in a general modulation scheme the Euclidean 
distances between various constellation point pairs differ and an error term ej depends on the symbol 
being transmitted. As a result, using ©, ©, and ©, the observation vector at node is 

y = H /i(uG © e) T + w. (8) 

Thus node has access to the likelihood p(y|u, e) and p(e) = YYj=iP( e j)> assuming the errors are 
independent. As shown in [3|, in order to avoid the propagation of errors occurring at intermediate nodes, 
node has to utilize the reliability information p(e). Then, the a posteriori probability of the source bit 
of interest, say u%, can be calculated by using the Bayes' rule: 

n 

p(« 1 |y) = a Yl P(yl u > e )II p ( e j)' ( 9) 

U2,—,Uk ei,...,e n j=l 

where a is a normalizing constant that does not depend on mj. The MAP estimate of u% at node is 
denoted by U\ and obtained as 

n 

ix\ = argmaxp(-u 1 |y) = argmax p(y|u, e) JJp(e 3 ), (10) 

141 Ul . , 

U2,— ei,...,e n 3=1 

which is the individually optimum detector for U\. As a result, for the optimal detection of u\, node 
requires the intermediate node reliability information vector: p e = [p ei . . . p £n ], where p e , = P(ej = 1) 
depends on the PMF of ej. We observe the performance of this detection rule in Section IV-A1 

The problem related to the MAP-based detection rule of (flOl) is the number of required operations, 
which grows exponentially both in the number of nodes k and the number of possible error events n. 
This is addressed in Section [TV] where we suggest a practical network decoding technique. 

C. Selective Network Coding 

The NC described in Section II-A is a static method in the sense that the generator matrix G is fixed. 
In static NC, node Vj always combines (network encodes) the symbols of a pre-determined set of users, 
even when it knows that the reliability for one of those users is low. When a symbol estimate with low 
reliability is combined with a symbol with high reliability, the reliability of the resulting network coded 
symbol is low. Thus, it is intuitive to expect some gains in performance by forcing the intermediate nodes 
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not to combine the symbols that have very low instantaneous reliability. In 11241 and 01, various forms 
of channel state information are used to determine thresholds for relaying decisions. In [17] and Il25ll , for 
relays assuming DF operation, successful decoding of channel code for a source is the required condition 
for combining its data in NC. Here, we propose a method called Selective Network Coding (SNC) that 
imposes a threshold on the reliability of the candidate symbols to be encoded at intermediate nodes that 
adapt DetF. In this way, any symbol that is sufficiently reliable is included in network encoding and the 
resulting encoding vector is appended to the transmitted packet so that node still has the instantaneous 
generator matrix G at the end of n transmissions. 

Let us demonstrate the operation under SNC on the sample network given in Fig. Q] For the first two 
time slots SNC is equivalent to NC since no combining of other nodes' symbols is the case. However, 
in the third slot, node 3 checks the reliability of the detection for U\ carried out following the first slot. 
Let us say it has observed an instantaneous SNR value of 7i_>3, which yields a probability of error 
equal to p ez = Q(^2'ji^ 3 ) in detection of u\ for BPSK modulation. Here, the Q-function is defined 
as Q(x) = -j= exp{— y)dz and the random variable 71^3 is exponentially distributed with mean 
value equal to average SNR 7 for a Rayleigh fading channel. The instantaneous error probability p ez is 
a measure of reliability for u-\_. This instantaneous p e3 value is averaged over 7^3 to set the threshold: 



where E 7l ^ 3 {.} denotes the expectation operator over 71^3. Therefore, node 3 uses the threshold value 
p t h 3 to check whether the detection at that instant is reliable. If p ez < p t h 3 , the detection is decided to be 
reliable enough and the combination u\ © w 3 is formed just in the way declared by the generator matrix 
G. Otherwise, node 3 modulates and transmits only its own symbol w 3 and appends this information to 
the corresponding packet. Similarly, in the last slot, node 2 checks the reliability of its own detection 
of u\ and forms either U\ © u 2 or simply transmits u 2 . Here, the reliability of Ui © u 2 is equal to the 
reliability of u-y. In general, there may be more than one symbol that an intermediate node should detect 
and combine according to G. In such cases the combined instantaneous reliability of an network encoded 
symbol at the time slot j can be easily obtained by 




(11) 



Pe 3 = 



(12) 



2 
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where A, denotes the set of sources for which vj should carry out the network coding, i.e., A, is the set 
of indices corresponding to the non-zero elements of the jth column of G, gj. The term Pj(ui ^ ih) in 
(fl"2l) is used to show the probability of error for detection of Ui by node Vj. 

Clearly, SNC inherently includes usage of adaptive generator matrices. The utilized generator matrix 
may assume in average a form dictated by some predetermined (and optimal if possible) linear block code 
structure like the ones that are to be discussed in Section Ull-B I Note that we do not claim the optimality 
of the proposed threshold. However, the bit error rate (BER) performance improvements are observed in 
Section IV-Dl 

III. Linear Block Codes Utilized as Network Codes 
A. Using Separation Vector as a Performance Metric 

Our goal is now to explore the error performance metrics for network coding/decoding described in 
Section HB Our basic figure of merit will be the diversity order corresponding to the source bit Uj, which 
is an asymptotic term defined for SNR tending to infinity: 

di = — hm ^—^ (13) 

SNR^oo log SNR 

giving information on the slope of decrease in logarithm of BER for «j, i.e., P^.^ u . for high SNR 
values. For conventional block coding, the average error performance over all data symbols is of interest. 
Therefore, for a linear block code whose coded symbols are transmitted over independent channels, the 
metric utilized for comparison is the minimum distance, which is equal to the diversity order 11261 . On 
the other hand, there is a vector of distinct minimum distances, i.e, separation vector, for data symbols, 
whenever we are interested in performance of individual symbols that originate from different source nodes. 
According to the results presented for suboptimal decoders in 030, 11201 , the diversity orders for symbols 
in some sample NC systems are still equal to the minimum distances in the corresponding separation 
vector in spite of the inherent error propagation problem. Using the soft decoding that we propose in 
(flOl ) and also authors analyze in lfl"2l . one should expect better performance and consequently diversity 
orders being equal to the minimum distances. A similar result is also shown for a simpler cooperative 
network with possible relay errors and the use of equivalent channel defined as the combination of the 
source-to-relay and the relay-to-destination channels Q. In Q, even a suboptimal detection rule utilizing 
this equivalent channel approach is shown to attain the achievable diversity order. As a result, supported 



12 

with intermediate node reliability information, the optimal rule of (TTOb given in Section III-BI also satisfies 
the diversity orders dictated by the separation vector whose entries are obtained according to Algorithm \T\ 
It should be also noted that since diversity order is an asymptotic quantity, the exact form of v is irrelevant 
to the procedure used for obtaining a diversity order value. On the other hand, it is wiser that each column 
gj of G is used as the encoding function for a Vj such that the jth entry is non-zero, Sj{ v j) 0- Otherwise 
possibly an extra relaying error is also included in the encoded data symbol. Therefore, v clearly affects 
the coding gain corresponding to the BER versus SNR curve of Ui. 

B. An Example of Close-to-Optimal Linear Block Codes: Greedy Codes 

In this study, we make use of some well-known linear block codes while constructing network codes 
that are to be used for the analysis of data rate and diversity orders for distinct symbols in Section IIII-CI 
and simulation of BER in Section |V] However, the cooperative network coded operation described in this 
work and the resulting performance figures for a unicast pair are more general and applicable to any linear 
block code like the maximum distance separable (MDS) codes detailed in the context of NC in ifTTl . 

In comparison with the network coded operation, we consider a case with no distributed coding (no 
network coding) among the nodes. For this no network coding scenario, we should also consider that our 
system model does not allow feedback of CSI within the network and that the average SNR values between 
all nodes are equal. If one intends to achieve higher diversity orders, two resources are available in such 
a scenario: (i) the temporal diversity resources over the faded blocks, (ii) the spatial diversity resources 
over the intermediate nodes. Here, it is seen that the source nodes must simply repeat their data instead 
of choosing a relay to convey their data which may possibly inject errors leading to worse performance 
than repetition. In conclusion, we call this method as the repetition coding scheme which is in fact a 
degenerate NC scheme with no cooperation hence with reduced spatial diversity resources. Following 
n transmissions, node combines the data received for each source symbol optimally to generate the 
detection results. 

On the other hand, with NC, we take the family of block codes known as greedy codes as an example. 
These (n, k, d) codes are selected with the following parameters: blocklength (number of transmission 
slots) n, dimension (number of unicast pairs) k, and minimum distance (minimum diversity order) d. 
Greedy codes are known to satisfy or be very close to the optimal dimensions for all blocklength-minimum 
distance pairs 11271 and can be generalized to non-binary fields [[281 for achieving higher diversity orders 
with NC as discussed in [fTvTl . Moreover, they are readily available for all dimensions (number of nodes) 
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and minimum distances unlike some other optimal codes. Hence, even in an ad hoc wireless network with 
time-varying size, any desired diversity order can be flexibly satisfied by simply broadcasting the new 
greedy code generator matrix G to be utilized in subsequent rounds of communication. 

As an example, let us consider a network that consists of k = 3 nodes transmitting their data symbols 
over GF(2). If a round of communication is composed of n = 6 transmission slots, we deal with codes 
of type (6, 3, d), which have a code rate of |. Starting with the generator matrix and scheduling vector 
corresponding to the repetition coding, we have 



G 



10 10 
10 10 
1 1 



v =[1 2 3 1 2 3] 



(14) 



It is easily observed that, since each data bit is transmitted twice over independent channels, this method 
satisfies only a diversity order of 2 for all bits U\, u 2 , and u 3 . In contrast, a diversity order of 3 for all 
sources can be achieved using NC, with the same code rate. As an example, the NC that achieves this 
performance can be obtained using the (6, 3, 3) greedy code, as follows: 



G 1 



10 110 
10 11 
110 1 



Vl = [1 2 3 1 2 3] 



(15) 



Clearly, without NC, the diversity order of 3 for all sources can only be achieved with rate ~. It should also 
be noted that greedy codes accommodate each unicast pair with equal diversity order due to the greedy 
algorithm utilized in their construction. Moreover, contrary to the findings in ifTTl . it is easy to obtain any 
required diversity order for any data bit even by using GF(2). The limitation is not due to the number of 
unicast pairs but due to the number of transmission slots in general. By increasing n, one can arrange and 
improve the diversity orders, if the transmissions to each node are realized over independent channels, 
which is a natural assumption for many wireless communication scenarios. If we need an increase in data 
rate, through a trade-off mechanism, we can assign decreased diversity orders to the lower-priority unicast 
pairs. This may be accomplished by omitting some columns of a greedy code generator matrix in order 
to decrease number of transmissions. The columns to be excluded can be decided by running Algorithm 
[T]in Section Hl-AI on candidate punctured generator matrices. As an example, the following punctured (5, 
3, 2) code is obtained by omitting the last column of Gi and has a data rate | that is higher than those 
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of above two codes: 



10 11 
10 1 
110 



v 2 = [1 2 3 1 2] 



(16) 



This punctured network code satisfies a diversity order of 3 for u\ and an order of 2 for both u 2 and u 3 . 
If Mi is of higher priority, this unequal error protection would be preferable especially when the higher 
rate of the code is considered. In case of larger diversity order requirement, d = 4 as an example, we 
may simply utilize the (7, 3, 4) greedy code with rate |. 



G, 



10 110 1 
10 111 
110 11 



,v 3 = [1 2 3 1 2 3 1] 



(17) 



A final problem is the selection of vector v. Our basic assumption is that v satisfies causality so that no 
intermediate node vj tries to transmit another node's symbol before hearing at least one copy of it. This 
causality problem can be solved trivially by using only systematic generator matrices. For the transmitting 
nodes corresponding to the non-systematic part of G, as described in Section IIII-Al one can select each 
entry Vj such that gj(fj) ^ for each column gj. For the columns that have more than one non-zero entry, 
a random selection between candidate v/s will merely affect the coding gains assigned to these nodes. 
As a result, one can force the number of transmissions of each node within a round to be equalized as 
much as possible for similar coding gain improvements of nodes. In the way exemplified in this section, 
one can choose a network code satisfying desired error protection properties for a determined network 
size with adequate data rate quite flexibly. 



C. Theoretical Gains in Rate and Diversity for NC 

In this section, we investigate the rate and diversity (asymptotic) gains of NC through use of the family 
of greedy network codes detailed in Section IIII-BL although the results are still valid for any other family 
of optimal or close-to-optimal codes. The availability of a greedy code for a given (k, d) pair is checked 
using [|29l . Fig. |2] shows the diversity gains attainable using greedy NC (with punctured codes in case no 
corresponding greedy code exists) with respect to the repetition coding scenario. The rate-diversity trade- 
off curves of both cases are plotted for a network of k = 3 nodes with increasing number of transmissions 
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and hence decreasing rate. We are interested in three types of network diversity orders; average, minimum 
and maximum, since the orders corresponding to each one of the three nodes may be unequal in general. 
The curves with no markers represent the (average) network diversity orders for both scenarios, which is 
defined as the arithmetic mean of orders for three nodes. For a rate of 0.43 bits/transmission, with greedy 
code (7,3,4), the network diversity order for NC is 4. The minimum, maximum, and average diversity 
orders are equal for this case. In contrast, the repetition scheme results in an average order of nearly 2.33 
with the worst node observing a minimum order of 2 and the best node a maximum order of 3, which 
would mean a high SNR loss asymptotically for all three nodes in the network. 

In Fig. [31 we now fix the desired network diversity order to d = 3 and observe the rate advantage of 
the NC for increasing network size. Note that for all cases diversity orders for k users are equal to 3. 
For a network of k = 25 nodes, the rate with NC is || (with greedy code (30, 25, 3)) and the rate of the 
repetition scheme is || (always equal to ~ for a diversity order of 3). The rate advantage ratio is then 2.5. 
In the asymptotic case, as k — > oo and hence as n — > oo, NC using optimal codes in construction will 
have a rate advantage converging to 3 since the rate for network coded case can be shown to tend to 1 
using the Gilbert- Varshamov bound ll30l for arbitrarily large n. In general, the rate advantage of NC over 
the repetition scenario becomes simply d, the desired network diversity order. As a result, increasing the 
network size improves the network coded system's efficiency in comparison to the repetition coding. 

IV. Sum-Product Network Decoder 

It is clear that the complexity of the optimal rule for decoding of any unicast transmission symbol 
Ui grows exponentially, since the number of additions and multiplications in (flOl) increase exponentially 
in the number of users k and transmissions n. Therefore, this rule becomes quickly inapplicable even 
for moderate-size networks. Recently the SP iterative decoding, which is often utilized for decoding of 
low-density parity-check (LDPC) codes, is suggested for decoding general linear block codes as well [|311 . 

Here, under the Rayleigh fading scenario detailed in Section III-BL we make use of SP decoding and 
compare its performance with that of the optimal one. The aim of the decoding operation is to produce a 
posteriori probabilities (APPs) for source symbols u\, . . . , To that end, we form a combined codeword 
[u\ . . . Uk c\ . . . c n ] and consider the parity check matrix for this codeword, which describes the underlying 
linear block code structure of the network code. On the Tanner graph, we add a variable node for each 
source symbol u^, i = 1, . . . , k and each coded symbol cj, j = 1, . . . , n. Afterwards, we add the check 
nodes which reflect the connections between the source and the coded symbols in the way described by 
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the parity check matrix. For the NC system given in (QQ), we refer to the graph presented in Fig. |4]for SP 
decoding at node 0. The parity check matrix for this system becomes: 
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u 3 
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(18) 



where G T denotes the transpose of G and I„ is the n x n identity matrix. For a regular LDPC decoder, all 
of the variable nodes are observed through the channel and corresponding to each channel observation an 
LLR is computed. For our case, the variable nodes u\, w 2 , and u 3 are not observed so the corresponding 
LLRs are set to 0. The channel LLRs for the remaining nodes (c l5 c 2 , c 3 , and c 4 ) cannot be calculated 
as in a regular LDPC decoder either, due to the intermediate node error events. Taking these errors into 
account by using © and ©, the channel LLR of c, is: 

p(Vj\cj = °) 



LLR c, 



In 



In 



In 



P(Vj\ C J 



(l -Pe 3 )p{yj 


Cj = 0) + p ej p{Vj 


Cj = 1) 


(l-p e .)p( yj 


&j = 1) + p ej p(yj 


Cj = 0) 



exp(LLR(e i )) exp(LLR(c i )) + 1 
exp(LLR(e i )) + exp(LLR(c i )) 



where 



LLRfe, 



In 



-^i and LLR(c,0 = In 
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P(Vj\cj 


= 0) 


P{Vj\cj 


= 1) 



(19) 



(20) 



where h* is the conjugated gain of the channel over which the modulated symbol Sj = fJ,(cj) is transmitted 
by node Vj and we use the fact that wj is Gaussian distributed (see Section Hl-BI) in obtaining LLR(cj). 
Given the channel LLRs, the SP decoder carries on iterations over the Tanner graph to generate the 
estimated LLRs for the source bits. If the number of iterations is fixed, the SP decoder utilized is known 
to have a complexity order of 0(n). In contrast, the optimal decoder has a computational load in the 
order of 0(2 n ), which makes the SP network decoder a strong alternative for increasing network size and 
number of transmissions. One may also note that the proposed scheme works directly with GF(g), q > 2, 
and constellations other than BPSK. The use of higher order fields and constellations would tremendously 
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increase the complexity of the optimal algorithm and make it impractical, whereas the SP algorithm would 
still operate with reasonable complexity. The number of iterations and other operational parameters for 
the SP decoder are given in Section IV-Cl where we show that performance figures close to that of the 
optimal one are possible for the network codes investigated herein. 

V. Numerical Results 

A. Sample Network-I: Simulation Results 

The results in this subsection are based on Sample Network-I of (OQ), consisting of only 4 nodes in 
order to observe the fundamental issues. For BER results, at least 100 bit errors for each data bit U\, 
1*2, and are collected through Monte Carlo simulations for each SNR value. In each run, data bits, 
intermediate node errors and complex channel gains are randomly generated with their corresponding 
probability distributions. The solid lines in Fig. [5] show the BER values for the optimal detector operating 
under the realistic scenario of intermediate node errors, whereas the dashed lines depict the performance 
of the genie-aided no-intermediate-error network with the same optimal detection. Finally, the dotted lines 
are for the detector that totally neglects possible intermediate errors. 

It is observed in Fig. [5] that different diversity orders for bits of different nodes are apparent for optimal 
detection under intermediate errors. The diversity order for U\ is observed to be 2 according to the slope 
of the corresponding BER curve. This is in agreement with the analytical results in Section III-AI where 
it was shown that an error event corresponds to at least 2 bit errors for the detection of u\ and u 2 . It is 
seen in Fig. [5] that the intermediate node errors cause no loss of diversity for u\ and u 2 , but an SNR loss 
of 1.5 dB. Hence the optimal detection rule of (flOl) is said to avoid the problem of error propagation in 
terms of the diversity orders. The loss for u 3 , whose diversity order is 1, with respect to the hypothetical 
no-intermediate-error network is around 2.5 dB. The performance deteriorates significantly for especially 
Mi and U2 when intermediate errors are neglected in detection (dotted lines), i.e., p es = p e4 = is assumed. 
Not only an SNR loss is endured but also the diversity gains for them disappear. 

B. Sample Network-II: Simulation Results 

Next, we verify the analytical results concerning the diversity orders for a set of three nodes operating 
under three different network codes constructed in Section IIII-BI Moreover, the unequal error protection 
performance of one of these codes is identified together with the rate advantage it provides. 
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The repetition method is represented by G and v in (fl~4l) . To construct Code-1 and Code-2, we make 
use of the greedy code of (fl~5T ) and the punctured greedy code of ([TBI respectively. Fig. [6] exhibits the BER 
curves for the repetition scenario with n = 6 transmissions (dashed lines), for NC scenarios with Code-1 
with n = 6 (solid lines) and Code-2 with n = 5 (dotted lines). The optimal detector of (flOl) is utilized 
for this simulation. Clearly, Code-1 has superior performance with an average network diversity order of 
3. However, the lower rate of Code-1 (and also repetition coding) in comparison to Code-2 should also 
be noted. For Code-2, on the other hand, bits u 2 and u 3 observe a diversity order of 2 while U\ observes 
an order of 3. With this unequal protection in mind, the average network diversity order for Code-2 is 
2+ g +3 ~ 2.33, which is higher than that of the repetition coding with order 2. In addition to improved 
diversity, Code-2 has also the advantage of increased overall rate and decreased decoding delay due to 
usage of 5 slots instead of 6. It is preferable especially for a network that puts higher priority on u\. 

C. Performance of the Sum-Product Decoding for Network Coded Systems 

In this section the performance figures for the SP iterative network decoder described in Section [IV] 
are presented in comparison with the optimal detection rule of (flOl) . which has an exponential complexity 
order. The network coded communication system of interest is given in (fl~5l) . The number of iterations for 
the SP type decoder is limited to 4 and no early termination is done over parity checks. Here, a minimum 
of 150 bit errors are collected for each data bit. 

In Fig. [7J we identify the fact that the SP decoder maintains almost the same BER performance as 
the optimal decoding rule. The SNR loss due to usage of the SP decoder is less than 0.1 dB for a BER 
value of 10~ 3 for all data bits. Achieving full-diversity with a linear complexity order, SP type decoding 
may serve as an ideal method for decoding for the network coded system of (fl~5l) despite the fact that the 
corresponding Tanner graph contains cycles. Results demonstrating the good performance of SP decoding 
were also reported previously in 011 - 11331 for graphs with cycles. In fact, one may realize that the length 
of the shortest cycle in the corresponding graph is 6, hence the graph is said to have a girth of 6. In ll33l 
within the context of sparse intersymbol interference (ISI) channels, it is shown that for any graph with 
girth 6, the performance of the SP algorithm is practically optimal. On the other hand, one may identify 
that for the family of greedy codes for k = 3 users with blocklenghts larger than 6, the girth of the 
corresponding graph will always be 4. Fortunately, it is also given in ll33l that the method of stretching 
on girth-4 graphs yields modified girth-6 graphs on which the SP algorithm evaluates the APPs for the 
data symbols with negligible performance loss. Further details on the girth profile and degree distribution 
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optimization procedures (like in [[341 for a greedy search of LDPC codes and like in 11331 for root-check 
LDPC code design) are out of scope of this work. 

D. Performance of Selective Network Coding (SNC) 

The selective network encoding operation defined in Section III-CI is applied in this section on the 
Sample Network-II of Section IV-BI The performance improvement for the selective encoding over the 
static (using fixed G with no selection of symbols to be encoded) encoding method is again shown 
using the SP iterative decoder of Section |IVj The instantaneous intermediate node error probabilities are 
compared with average error probabilities (dictated by G) and data of the nodes whose error probabilities 
are below the corresponding average values (thresholds) are combined by the intermediate node. In Fig.[8l 
we observe that SNC offers an SNR improvement of 0.6 dB for BER set to 1CT 3 over the static NC method. 

E. Performance of NC under Slow-Fading Channel Model 

All the discussion and the results presented up until this section rely on the assumption that all channel 
gain coefficients related to the observations at node are independent. Hence a block fading model over 
time slots is utilized. However, it is also possible under many communication scenarios that the variation 
of a channel gain coefficient is not rapid enough for such an assumption. Then it is also possible that all 
transmissions from a selected source node to node observe the same fading condition leading to the 
degradation in BER performance due to loss in diversity. Therefore, we finalize the numerical results by 
providing the BER curves of NC and repetition coding under the assumption that within a round of n 
transmissions, only the transmissions from distinct source nodes observe independent fading, i.e., hj and 
h m are independent if vj ^ v m and otherwise hj = h m . In Fig. |9l we investigate the BER curves for the 
Sample Network-II operating under this slower fading assumption. The repetition coding is represented 
by (fl4l) and NC is realized by (fT3T ) . It is seen that the repetition coding merely results in a diversity order 
of 1 for each symbol as expected. On the other hand, NC yields an order of 2 via the cooperative diversity 
obtained due to intermediate nodes transmitting over independent channels. The SNR losses incurred by 
not utilizing NC are shown to further increase in great amounts for this slower fading channel scenario. 

VI. Conclusions 

We formulated a NC problem for cooperative unicast transmissions. A generator matrix G and a 
scheduling vector v are used to represent the linear combinations performed at intermediate nodes. We 
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presented a MAP-based decoding rule utilizing G, v, and the error probabilities at the intermediate nodes. 
A method for obtaining the performance determining parameter as the diversity order for individual source 
nodes is proposed for any given G over the corresponding separation vector. Through simulations we 
showed that our decoding rule, using reliability information for the network coded symbols, avoids the 
diversity order losses due to the error propagation effect. We presented design examples for network 
codes via greedy block codes, which may also provide unequal diversity orders to nodes with proper 
puncturing. Over given design examples, we obtained the rate-diversity trade-off curves and the rate 
advantage realized by using NC with respect to the no NC case. Moreover, the SP iterative network 
decoder with linear complexity order is proposed and shown to perform quite close to the optimal rule. 
Furthermore, the selective NC scheme combining only the reliably detected data at cooperating nodes is 
shown to yield additional coding gains. Identifying gains of NC for purely random G matrices in large 
networks, studying the effects of imperfect information on channel gains and relay error probabilities 
will be addressed in future work. Finally, it would be also interesting to operate suggested wireless NC 
methods under asymmetrical channel gains, which can be more realistic for ad hoc networks. 
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Fig. 1. Sample network coded transmission scenario 
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Fig. 2. Network diversity orders for greedy NC and repetition coding for given data rates. 
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Fig. 3. Increasing rate advantage of greedy NC for increasing network size. 



+ 




+ 




+ 




+ 




u 3 



C2 



C3 



c 4 



Fig. 4. Tanner graph for network coded system of (Q} 
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Fig. 6. BER performance for repetition coding and NC with greedy codes. 
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Fig. 7. BER curves for the individual MAP decoder of UOt and the SP iterative decoder 



2S 




10 1 1 1 1 1 1 

2 4 6 8 10 12 

E/N (dB) 



Fig. 8. Selective and static network encoding BER curves 
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Fig. 9. BER vs. SNR curves for slow fading channel 



